euler_028 = sum $ take (1+4*((1001-1) `div` 2)) $ scanl1 (+) $ 1 : (concat $ map (replicate 4) [2,4..])